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The transition probability for time-dependent unitary evolution is invariant under the reversal 
of protocols. In this article, we generalize the expression of microscopic reversibility to externally 
perturbed large quantum open systems. The time-dependent external perturbation acts on the 
subsystem during a transient duration, and subsequently the perturbation is switched off so that 
the total system would thermalize. We concern with the transition probability for the subsystem 
between the initial and final eigenstates of the subsystem. In the course of time evolution, the energy 
is irreversibly exchanged between the subsystem and reservoir. The time reversed probability is given 
by the reversal of the protocol and the initial ensemble. Microscopic reversibility equates the time 
forward and reversed probabilities, and therefore appears as a thermodynamic symmetry for open 
quantum systems. 

PACS numbers: 05.30.-d,05.70.Ln 



I. INTRODUCTION 

The time reversal invariance of the equations of mo- 
tion amounts to universal symmetry of fluctuation the- 
orems and related equalities [H-GJ] for the fluctuation 
of entropy production and particle current for meso- 
scopic systems, and thus plays fundamental role in the 
nonequilibrium statistical mechanics. The symmetries 
connect theprobabilities of positive and negative entropy 
production]^, [2(|, which are calculated from the time 
forward and reversed transition probabilities for quan- 
tum systems P. 15, Il3l. Il5l - fl7j . In the classical Markovian 
stochastic dynamics, the conditional probability func- 
tional of the trajectories satisfies a symmetry expressed 
by the probability functionals and heat. This relation is 
also called microscopic reversibility @, and generalized to 
the quantum open systems by concerning the heat calcu- 
lated from set of transitions of the reservoir states [lol [ill ]. 
On the contrary, we are interested in the transition prob- 
abilities between the system states instead of the statis- 
tics for dissipative quantities, and give an expression of 
microscopic reversibility. It is known that the induced 
absorption and emission by an external electric field are 
equally probable within the realm of the Fermi's Golden 
rule provided that the initial and final Fock states are 
exchanged. Absorption process can be seen as the time 
reversal of the corresponding emission process. 

In this article, we show that a symmetry similar to 
that of the induced absorption holds for the macroscopic 
time-dependent open systems where the system couples 
to the reservoir. Since the expression appears as a generic 
equation for the forward and reversed protocols, we call 
it as microscopic reversibility. 

This article is organized as follows. First, we describe 
our model and forcing protocol. And an expression of 
microscopic reversibility for open systems is derived in 



Eq. (ll4[) . Then the microscopic reversibility is numeri- 
cally verified. 



II. MODEL 

Let us consider a finite system interacting with a 
macroscopically large reservoir at an inverse tempera- 
ture p. The system is externally controlled by a time 
dependent parameter X(t), which is for example a spring 
constant for the case of a harmonic oscillator. Therefore 
the total energy change of the system is caused by the 
external work done and the heat flow from the reservoir. 
The total Hamiltonian is 



H{t) = H s (X(t))+H r + H s 



(1) 



where H s (X(t)), H r , and H sr are the Hamiltonians of the 
system, the reservoir, and the interaction between them, 
respectively. 

Let us prepare the initial state as 



P(0) 
P.(A(0)) 



p s (A(0)) ®p r ; 
e -0H a (\(O)) 



(2) 



which is the product of the canonical ensembles of the 
system p s (X(0)) and of the reservoir p r at the same in- 
verse temperature /3. The partition function of the sys- 
tem is 



Z(X(0)) = Tr s e 



(3) 
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Through out this paper, we assume that the interaction 
energy H sr is small compared to the energy of the sub- 
system H s and the reservoir H r . This assumption is rea- 
sonable for macroscopic systems, since H s and H r are 
the bulk energy, while H sr would be proportional to the 
surface area. Note that the assumption is different from 
the weak coupling limit. Namely, the weak coupling as- 
sumes that interaction is negligible, and the perturbative 



2 



analysis is available. On the other hand, in our case the 
interaction is not necessarily vanishing, while we require 
that the ratio between the interaction energy and bulk 
energy is negligible. 



III. TWICE MEASUREMENTS SCHEME 



describes the transition of the system state, and 
satisfies 



W)) 
\m(T))(m(T)\U\n(0)) = 1. 



(10) 



The twice measurements scheme consists of the initial 
and final observations, (i) and (ii). 

(i) At t = 0, we measure the energy operator H s (X(0)), 
and gain an eigenenergy E n (0). The system state 
becomes the corresponding eigenstate ".(0)}, 



tf s (A(0))|n(0)) = E n (0)\n(0)}. 



(4) 



Subsequently the total system unitarily evolves un- 
tiH = T as 

e -/9^»(o) 

p{T) = [/(|n(0)) (n(0)| ® Pr )U+. (5) 

Here U = T{e~^ fo dsH ( s )j is the unitary time evo- 
lution operator. For the time evolution, we require 
that the external forcing is switched off well in ad- 
vance t — T and the total density matrix would 
relax to an equilibrium state at t = T, 



Up{0)U+ S p.(A(T)) ® Pr ] 

e -PH s (\(T)) 



Ps(HT)) 



Z(X(T)) 



(6) 



Here H S (X(T)) is the corresponding energy opera- 
tor at t = T and Z(A(T)) is the partition function. 
Eq. comes from the equation at the level of ma- 
trix elements [2 5] 



Tr r Up(0)U + = p s (X(T)) 



(7) 



and the assumption of smallness of the interaction 
energy, and holds when acting on the local state of 
the subsystem as in Eq. ([T4|) . 



(ii) At t = T, we measure the energy H S (X(T)) and ob- 
tain some eigenenergy E m (T) with the correspond- 
ing eigenvector \m(T)). Regarding the reservoir, we 
don't perform any measurements. The transition 
probability that the initial and final system states 
are |n(0)) and |m(T)) is then 



P F (\n(0)) -> \m(T)}) 

e -/9B„(0) 



Tr r {(m(T)|C/|n(0)) 



Z(A(0)) 



p r (n(0)\U+\m(T))}. 



Here the Kraus operator 

1 PE n (0) 
A = p 2 

Vz(xW) 



(8) 



m{T)\U\n(0)) (9) 



Note that the matrix elements such as 
(m(T)\U\n(0)) contain the reservoir variables. 
The second equality follows from the complete- 
ness J2 m l TO (^))( TO (^) = 1) an( i the unitarity 
U+U = 1. 

Similarly the probability of the time-reversed dy- 
namics is calculated as well. Firstly, let us define 
the reversed dynamics by reversing the time de- 
pendence of the system Hamiltonian in Eq.(2) as 
H S (X(T — t)), i.e. the reversal at time T. Also, we 



start with the initial state 0- 



€5 p r & 



z(\(T)) 

where & is the anti unitary time reversal opera- 
tor of the total system. The partition function is 
defined as 



Z{X{T)) = Tr s e 



-0H s (\(T)) 



(11) 



(i-2) We measure the energy 67? S (A(T))9^ 1 at t = 
and consider the case that the eigenenergy E m {T) 
corresponding to (i), and the system state becomes 
0|m(T)}. 0|m(T)) is proportional to |m(T)} when 
the energy E m (T) does not degenerate. 

(ii-2) At the final time t = T, we again measure the 
energy 07J S (A(O))0 -1 and obtain the eigenenergy 
E n (0). The system state is 6|n(0)). 

Then the probability that the initial and final states 
are 0|m(T)) and 0|n(O)) for the reversed dynamics 
is 

P R (@\m(T)) -> 0|n(O)» 

, „ — f3E m (T) , 

= TT I {{n(0)\eU + e\m(T))——— Pr {m(T)\eue\n(0))}. 



Z(X{T)) 



(12) 



Here (n(0)| &U+Q\m(T)) = (m(T)\U\n(0)) is the 
inner product of 0|n(O)) and U + Q\m(T)). 



IV. MICROSCOPIC REVERSIBILITY 

Now let us show the relation between probabilities of 
the forward and reversed transitions 

Pf (MO)) -»■ |n(T)» = P R (Q\n(T)) -> 9|m(0)». (13) 

Here the difference between the forward and reversed 
probabilities goes to zero in the macroscopic limit. For 
this purpose, we apply the methodology developed in 
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Refs.fiol. HTj and thermalization property Eq.©. We 
rewrite the forward probability as 



Tr r p|(n(0)|- 



|m(T)» 
1 

— i — i 

Pr 2 PrPr 2 (m(T)\U 



e-^ H '^U+\m{T)) 



■i/3ff s (A(0)) 



n(0))/9= 



= Tr r (n(0)|C/ + ({[/p s (A(0))v?[/+} (Os (A(T))-V-2 ) 
p s (A(T))i|m(T))p r (m(T)|p s (A(T))i 
( Ps (A(T))-v;^{C/ Ps (A(0))^|(7 + })C/|n(0)) 
Tr r (n(0)|£/ + 7=^= ^^ A ( T >V(T)y r 



VzWT) 



<T)\ 



1 



1 



Z(A(T)) 
P H (0|m(T)) 



e-^ m(T) Tr r (n(0)|C/ + |m(T)}p r (m(T)|;7|n(0)) 
8|n(0)». (14) 



The first equality derives from the cyclic property of the 
trace. In the second equality, U + U = 1 is inserted. Also, 
with the use of the relaxation property of the density 
matrix Eq.tfTBf which comes from Eq.0 and the small- 
ness of interaction energy H sr with respective to the 
bulk energy, the contents of the curly brackets are just 
the inverse of the remaining quantities in the brackets 



V. NUMERICAL DEMONSTRATION OF 
MICROSCOPIC REVERSIBILITY 

In this section, we numerically show the microscopic 
reversibility Eq. (IT4"|) ■ See also the detailed numerical ver- 
ification in Ref. [221] for various parameters and measure- 
ment basis. We consider N = 8 site spin chain in the 
spatially inhomogeneous time dependent magnetic field. 
Note that thermalization property has been observed in 
relatively small system sizes [23J, |24|. In the context of 
quantum derivation of thermal state, the dimension of 
the Hilbert space, which exponentially depends on the 
system size plays essential role[25|. The Hamiltonian is 
given as 



H(t) = H s {t) + H r {t) + H a , 



N.-l 



N.-l 



N s — 1 



H s (t) = -J s ]T <r|<r| +1 + a s £ of + h(t) °i 



N 



N 



N 



H r = -J r J2 °i°i+l+<*r J2 "i+Tr E °i 
j=N 3 + l j=N 3 + l j=N 3 + l 

H sr = -Jsr&N 3 a N a +l + a sr&N s + lsr^N 3 , (15) 

where crj is the i component of the Pauli matrix. The 
sites j = 1,2(N S = 2) and 3 < j < N are regarded 



P F (M0)>->|n(t)>) 
0.010 



0.008 



0.006 



0.004 



0.002 - 




FIG. 1: The time dependence of the forward and re- 
versed probabilities Pp(|m(0)) — » |n(t)))(blue line) and 
Pii(Q\n(t)) — > O|m(0)}) (black line) is shown where the ini- 
tial and states are the eigenstates with the first to second 
largest eigenvalues. 



as the subsystem and reservoir. The external magnetic 
field h(t) = f - f tanh(/i(i - r)) satisfies h(0) = B and 
switched off after t 



t. The initial state is prepared as 



p(0)=p s (0)pr 

1 -BH r 
Pr = 7J-e H . 



(16) 



We calculated the forward and reversed transition prob- 
abilities in Fig.l. The parameters are chosen as follows. 
The exchange interactions at each site are J s = J r = 
J sr = 1, the x component of the magnetic field is ex- 
pressed by the parameters a s = a r — 1, a sr — 0.2, and 
similarly the z component of the magnetic field is given 
by h(t) with B = 3, fj, = 5, 7 r = 1, and "f sr — 0.2. The 
inverse temperature is /3 = 0.01 and the switching time 
is t = 5. 

The unitary time evolution U = T{e~* fo H(t)dtj 
is discretized as e -i*tH(NAt) e -iAW((N-i)At) . . . 

e -iAtH(At) e -iAtH(0) with the time gtep At = o 05 which 

is much shorter than the time scale of external perturba- 
tion. Initially, we measure the eigenstate of H s (0) cor- 
responding to the largest eigenvalue. The measurement 
basis at t — T is chosen as the eigenstate of H s (T) which 
corresponds to the second largest eigenvalue. Note that 
there are no crossing of time dependent energy levels of 
the subsystem, and Eq. ([TB]) holds also for the transitions 
to other three eigenstates. It is remarkable that the time 
dependence of the forward and reversed transition prob- 
abilities are similar as a function of time even after the 
quench at t = r. 
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VI. SUMMARY 

We have derived an expression of microscopic re- 
versibility for macroscopic quantum open systems. The 
microscopic reversibility is trivial for the transition prob- 
ability of unitary time evolution for the total system. 
Similar reversibility is well-known for the transition prob- 
abilities of induced absorption and emission under the 
influence of an electric field in equilibrium. The micro- 
scopic reversibility is regarded as a generalization to the 
case of generic time-dependent perturbation. Therefore, 
it is a symmetry holds in generic macroscopic quantum 
open systems. The microscopic reversibility is numeri- 
cally verified for a spin chain with time-dependent per- 
turbation. In the context of the quantum generalization 
of the reversibility for the classical conditional probabil- 
ity functional, Ref.(l0| derives another symmetry for the 



quantum trajectories. The main difference from Ref.[10( 
is the quantity which we measure, number of measure- 
ments, and the definition of the reversed process. Here we 
pursue and measure the system states. This is a useful 
property of the present scheme, since the system vari- 
ables are expected to be much easier to measure with 
sufficient accuracy compared to those of the spatially ex- 
tended large reservoir. 
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